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Abstract 

The dynamical symmetries of the Kratzer-type molecular potentials (general- 
ized Kratzer molecular potentials) are studied by using the factorization method. 
The creation and annihilation (ladder) operators for the radial eigenfunctions sat- 
isfying quantum dynamical algebra SU (1, 1) are established. Factorization method 
is a very simple method of calculating the matrix elements from these ladder op- 
erators. The matrix elements of different functions of r, r-3~, their sum T\ and 
difference T2 are evaluated in a closed form. The exact bound state energy eigen- 
values E Ut £ and matrix elements of r, r-4-, their sum T\ and difference T2 are 
calculated for various values of n and t quantum numbers for CO and NO di- 
atomic molecules for the two potentials. The results obtained are in very good 
agreement with those obtained by other methods. 
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1 Introduction 



Exactly solvable models in Physics have been generating a renewed interest, because 
of the fact that they can be solved in terms of creation and annihilation operators 
by means of factorization method [1 - 5]. In the year 2002, this present method was 
proposed in the two related works: the ladder operators for the modified Poschl- Teller 
potential and the Morse potential were obtained [3J. With factorization method, the 
ladder operators of a quantum mechanical system with some important potentials like 
Morse potential, Poschl- Teller potential, the pseudoharmonic potential, the infinitely 
square- well potential and other quantum systems have been established [3j H]. For 
comprehensive review (see Dong and references therein [2]). 

Algebraic methods underlying Lie symmetry and its associated algebra have been 
widely used to study many of these exactly solvable potentials. In 2007, Rasinariu et 
al. [5] gave a review of the progress made so far in solving exactly solvable problems in 
quantum mechanics, by connecting supersymmetry and spectrum generating algebras 
through the property of shape invariance. The solutions of exactly solvable models can 
be achieved by dynamical algebraic approaches [1-11]. In 1991, De Lange and Raab [6] 
presented operator methods with shift operators (that is, raising and lowering operators) 
for the Hamiltonian of exactly solvable models. 

For the factorization method, we are going to adopt Dong's approach [2] by finding 
the ladder operators with the following properties: 

K±$ n (x) = k n <S> n+1 (x). (1) 

In this case, we are seeking for the ladder operators of the form 

K£ = A$(x)± + I%(x), (2) 

these ladder operators depend on the physical variable x, which is different for different 
quantum systems. With equation(JTJ), the following expressions are obtained: 

K n+i K t®n{x) = k+K- +1 ^ n (x), K+K- +1 ^ n+l (x) = k+K- +1 $ n+l (x), (3) 

from which we find that the products of the operators K~ +1 K+ and K^ +1 K~ acting on 
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the given wavefunctions $> n (x) and $ n+ i(x), respectively, have the same constant 

c n = k+K~ +v (4) 

The SU(1, 1) algebra has useful applications in Physics [21 [Til EH E]. The algebra 
of the group SU(1, 1) is used to generate the energy spectra while the representation 
matrices of the group could be used to calculate time dependent excitations of the bound 
states and the scattering states respectively 

For diatomic molecules, the generalised Kratzer potentials is considered as one of the 
molecular potentials. These potentials have the general features of the true interaction 
energy, interatomic, inter-molecular and dynamical properties, their wavefunctions do 
vanish at the origin like Morse potential. This added advantage make these potentials 
important in Molecular Physics, Chemical Physics, Solid State Physics. These potentials 
(Kratzer and modified Kratzer molecular potentials) that is being considered in the 
present work are important molecular potentials which describe the interaction between 
two atoms . These potentials have been used extensively to describe the molecular 
structure and interactions and have been receiving much attention in the history of 
quantum chemistry for some decades ago [121 EE I2Q1 (2TJ [23] . 

It is the purpose of this paper to study the dynamical symmetries of the Kratzer-type 
molecular potentials (generalized Kratzer molecular potentials) and to establish that the 
generators obtained are an S77(l,l) dynamical group, and to obtain some numerical 
results for the energy eigenvalues and matrix elements for some diatomic molecules. 

The paper is organized as follows. In Section 2, we study the exact solutions of the 
Kratzer-type molecular potentials. Section 3 contained the construction of the creation 
and annihilation operators, the commutation relations of some of these operators and the 
matrix elements of some related functions r and r4-. Also, in Section 4, the numerical 
calculations of the energy levels and matrix elements for some diatomic molecules for 
the Kratzer and modified Kratzer potentials are given. We conclude with Section 5. 
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2 Exact solutions of the Kratzer-type molecular po- 
tentials 



Consider the molecular potential (the generalized Kratzer potential) of the form [TTJ [T8| 

using 

. . a b 

V(r) = - + - + c. 5 

This potential is of great interest because, it is a generalised form of the molecular 
potential which can be used to generate other forms of the Kratzer-type molecular 
potentials: 

• Standard Morse or Kratzer-Fues potential 

nr) = -D, - f) (6) 

where Dq is the dissociation energy between two atoms in a solid and ro is the equilibrium 
intermolecular separation. On comparing with equation (5), we have a = —2D Q r ,b = 
D Tq, c = and this potential has minimum value to be —-Do at r = r [2, 12 - 14, 17 - 
23]. 

• Modified Kratzer potential 



where, a = — 2DqTq, b = DqTq and c = Dq. Modified Kratzer potential is Kratzer 
potential shifted by the amount D [13 [TBI US ED] ■ The graphs of the two potentials 
(equations 6 and 7) are shown below in Figure 1. 

Figure 1. Shapes of the Kratzer and modified Kratzer potentials for CO and NO 
diatomic molecules. 

Consider the motion of a particle in a spherically symmetric potential V(r) 

- ^A*(r, 9, <p) = [E — V(r)]*(r, 9, 0) (8) 

where the potential V(r) is taken as in equation (jSJ). We seek for the wavefunctions of 
the form 

*nAm(r, 0, 0) = RnAr)Y e m (0, 0) (9) 



which reduces equation (jHJ) into the radial and angular wave functions as: 



d 2 R nl (r) 2 dg^ (r) 



2/i 
F 



— 



a 6 > 

- + +C 



and 



L 2 F7( 



h 2 £ (£ + l)Y e m (9, 



£(£+1) 



0. 



^ (r) = (10) 
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where R n ^(r) and Y™(9, 4>) are the radial and angular solutions of equations (ITOland fTT]) 
respectively. For ^ n ^ m (r, 9, 0) to be finite everywhere, R nj g (r) must vanish at r = 0, 
that is, R n> £ (0) = 0, then R Ut g (r) is a real function. 

For the bound state energy eigenvalues for this quantum system, the following di- 
mensionless abbreviations are introduced : 

—/J, \ 2 a I a i\ 



2 -(£-e) 



r; o; 



a; 



+ £{£+!). (12) 



^h 2 (E-c)J ' ^^'^ ^ 
The above equations (TT2k ) and (|12b ) are chosen in order that the acceptable bounds state 
solutions are obtained, this can only be possible if E < (otherwise, continuum state 
solutions will be obtained). /3/ in equation ffl"2b ) represents the usual centrifugal term, 
that is, /3g (fit + 1) = ^ + £ (£+1), where £ is the angular momentum quantum number. 
These substitutions allow us to obtain the following hypergeometric-type equation. 

d 2 



dp 



■Rn,£ (p) + --7-RnJL (p) + 

p dp 



1 a &(& + !) 

a n 



Rn,t (P) = 0. 



(13) 



4 p p 

This differential equation has an irregular singularity as p — > oo, where its normalized 
solutions in bound states behave like exp (p). It further has a singularity at p — > 0, where 
Rn, i (p) ~ f 1 - Then, the ansatz for the wave functions which is a physically acceptable 
solution for R nt £ (p) can be expressed in the form 



Rn,l (p) = N n/ e-»' 2 p^G{ P ) 



(14) 



and therefore, equation fflBl becomes 
d 2 G(p) 



P- 



[(2/3/ 



p}^--[P e + l-a]G( P ) = 0. 



(15) 



dp 2 LV / rj dp 

This is the associated Laguerre differential equation (Kummer equation) [TU [23], the 
solution of equation f[T5"j) which is regular at origin (regular at r = or p = 0) is the 



degenerate hypergeometric function 

G{p) = ^(ft + 1 - a, 2/3, + 2; p). (16) 

For large values of p, this solution diverges as exp(p), thus preventing normalization, 
except for when 

{3g + 1 — a = —n r = —n; n r = 1, 2, . . . (17) 

becomes a polynomial. 

Therefore, the solution for the radial equation for this generalized Kratzer-type 
molecular potential is 

R n ,i (p) = N n/ e~ p/2 p^ iFa(-n, 2/3, + 2; p), (18) 

where the normalization iV nri , is determined from the requirement that 

poo 

/ I R n/ (p) | 2 r 2 dr = 1. (19) 
Jo 

By using the expression that relates the associated Laguerre functions with the con- 
fluent hypergeometric functions 

7 \m ! 



(7 + m) ! ' 



iFi(-7, m + 1; z) 

together with the following important formular [TH [23] 

(a + 2ra)r(a + n) 



/ e-*x a L a - l {x)L«-\x)dx 
Jo 



n ! 



and by substituting 7 = £ 2 , the normalized radial wave function is obtained as 



Rn,e (r) = N n/ e" V<2^ +1 (£r), 



where 



n! 



2 (n + ^ + l)r(n + 2/3, + 2) 
The corresponding eigenvalues are : 



E„ 



-pa 



2h 2 [n + /3, + 1] 



2 +c, 



where /3, is obtained as the positive root of equation (fl2"b ) given as 

1 



-i + v/(^ + i) 2 + ^ 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 
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3 Construction of the creation and annihilation op- 
erators 



The ladder operators can be generated directly from the eigenfunction with the factor- 
ization method as shown in Dong [2] and references therein. We shall find the diferential 
operators L± with the following property: 

L ± R n/ (r) =£ ± R n±l/ (r), (26) 

the operators of the form 

L± = A ± {r)^- + B ± {r) (27) 
dr 

which depend only on the physical variable r are to be obtained. 

The action of the differential operator j- on wave functions (1221) gives: 

^-R n/ (r) = (r) + OLr^ (r) + N n/ e~* -^ +1 (£r). (28) 

dr 2 r dr 

The expression above is used to construct the ladder operators L± by using the recur- 
rence relations of the associated Laguerre functions in order to find the relation between 
R n ,£ ( r ) an d R n +i,e ( r )- To find these, the following recurrence relations of the associated 
Laguerre functions are used [HJ 121]: 

,r—T a (^- I nL n( X ) ~(n + tt)^-l(^) (0Q\ 

dx n{ ] \{n+ l)L% +l (x) -(n + a+1- x)L%(x) 1 } 

and the creation and annihilation operators are obtained as: 

d <f~T d £t 

C- = -r— - Y + n + p e , £+ = r— - ^ + n + + 2, (30) 

where n is the number operator with the property 

nR n:£ (r) = nR n/ (r). (31) 

The action of the creation and annihilation operators on the radial wavefunctions 
R n ,e (r) gives the following properties 

C±R n/ (r) = £±Rn±i,i (r), (32) 



where 



\ 



n(n + p e )(n + 2fa + 1] 
(n + /3 e + 1) 



\ 



(w + l)(n + & + 2)(n + 2& + 2) 
+ + 



(33) 



On studying the dynamical group associated to the annihilation and creation oper- 
ators £_ and £ + and based on the results of equations (132]) and fl33l . we can evaluate 
the commutator [£_, £ + ] as: 



where 



and £ n is defined as 



[£_, £ + }R n/ (r) = 2i R n/ (r) 

e = (n + p t + 1) 

£o = (n + Pe + l). 



Thus, operators £ T and £$ satisfy the following commutation relations: 

[£o, £ T ]R n ,i (r) = T£ T R nT i,i (r). 

The action of £ + on the radial wavefunction R n ^ (r) gives 

(£ £+ - £+£>o) Rn,e (r) = £+R n ,e (r) 
= £ (£+R n/ (r)) - £ (£ + R n/ (r) 
= (£ + l)£ + R n/ (r). 



(34) 
(35) 
(36) 

(37) 



(3* 



If £+Rn, i if) is non-zero, then, it is an eigenfunction of £q with eigenvalue (£o + 1). 
Thus, the effect of £ + is to raise the eigenvalue by one unit. Similarly, the action of £_ 
on the radial wavefunction R n j (r) gives 



(39) 



£ £- - £-£ ) Rn,£ (r) = -£- R n ,e (r) 
= £ {£- R n ,i (r)) - £ (/L R n/ (r)] 
-- (£ - 1)£_ R n , e (r), 



and if £+R n , e (r) is non-zero, then, it is an eigenfunction of £ with eigenvalue (£ — 1) 
(this is the reason for calling £ + and £_ raising and lowering operators respectively). 
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For the Hermitian operators, we define the operators as follows: 
-(4 + 4), £ v = - 



4 = -(4 + 4), 4 = -(4-4), 4 = 4, (40) 



we obtained the following commutation relations 

The Casimir operator (25] can be expressed as 

CR n/ (r) = (4(4 - 1) - £+£-) Rni (r) 



= (4(4 + 1) -£-£+) Rn,e (r) (42) 
= Pt (ft + l)Rn,i (r). 

The Casimir operator C now satisfies 

4 4] = 4 4] = 4 4] = 4 4] = o, (43) 

the operators £±, £ x , £ y , £ z and 4 satisfy the commutation relations of the dynamical 
group SU(1, 1) algebra, which is isomorphic to an £0(2, 1) algebra (i. e. 577(1, 1) ~ 
50(2, 1). The commutation rules are valid for the infinitesimal operators of the non- 
compact group SU (1,1) [TOf [26]. 

These relations coincide with the formulas that define the action of the infinitesimal 
operators £± and 4 of the dynamical group SU (1, 1) on a basis | j, k) of the irreducible 
representation D + (j) belonging to the discrete positive series in an abstract Hilbert space 
[21 [TOl [HI [26]. The eigenvalues have the ground state and therefore, the representation 
of the dynamical group 577(1, 1) belongs to D + (j), using the Dirac notation | j, k) = 
Rni (r) : 

CJj,k)=j(j + l)\j,k) 

4 I j, k)=£o\ j, k) = k\ j, k) 



4 I j, k) 



k(k±lf-j(j + l)(k±l) 



(44) 

j,k±l) 



k 

k = -j + n, n = l,2, ...,j<0. 
Furthermore, the following expressions can be easily obtained from the operators 4 
and 4 as follows: 

r = i[24-(4 + 4)] 
d * 1 ~ ( 45 ) 

-H<a-^)-i. 



With these, the matrix elements for r and r j- are obtained as follows: 

(R m/ (r) | r | R n/ (r)) = - [(n + f3 e + l)8 m>n - £ + 5 m , n+1 - £_5 TOi „_i] (46) 

and 

(^ m ,< (r) I r— | (r)} = — 5 mri+1 - — 5 m , n _i - <5 m ,„. (47) 
From equations (JIB]) and (jlTl) . We can deduce the following relations: 



(r) | r | (r)> + (R m>i (r) | r— | (r)> 

1 3£_ 
= (n + A )5 m , n - -£ + 5 m ,„ +1 - — ^ 

and 

£(-RrM (0 | r | i? n ^ (r)) - (R rrh i (r) | r— | R n/ (r)) 



(46 



(49) 

these relations form a useful link for finding the matrix elements from ladder operators. 
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(n + A + i)l,„--i + W-YVi> 



4 Numerical calculations of the energy levels and 
matrix elements. 

4.1 Energy eigenvalues of the Kratzer and modified Kratzer 
potentials. 

The energy eigenvalues for the Kratzer and modified Kratzer potentials are obtained 
respectively, as : 



ft ra + + 1 



and 



where 



rpMK 2/iD Q r 

2ft n + + 1 



2 



(52) 



In this work, energy eigenvalues for CO and iVO diatomic molecules for the various 
values of n and £ are obtained by means of the factorization method (FM) (equations 
(50) and (51) ) with the parameters given in Table 1. The results obtained are compared 
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with other results obtained by using: AIM method |12j ; Nikiforov-Uvarov (NU) method 
[15] and Exact quantization rule (EQR) method [20J. 

Table 2 shows the exact bound state energy eigenvalues of the Kratzer potential 
for CO and NO diatomic molecules for various values of n and I using equation (50) 
obtained by factorization method (FM) and other results obtained by using AIM [12J and 
EQR [20] methods. Similarly, Table 3 shows the exact bound state energy eigenvalues 
of the modified Kratzer potential for CO and NO diatomic molecules for various values 
of n and I using equation (51) obtained by factorization method (FM) and other results 
obtained by using AIM [12] and NU [15] methods. 



4.2 Matrix elements of r and r jj- for the two potentials. 

The matrix elements for r and t4- are given as: 



(Rnj (r) | r | R n/ (r)) 



1 
1 



{n + p t + 1) 



(n + 2)(n + /3 £ + 3)(n + 2/3 £ + 3) 
(n + Pt + 2) 



\ 



(n-l)(n + ^-l)(n + 2^) 



{n + Pt ) 



(53) 



and 



(R n/ (r) | r^j- | R n i (r)) 



1 
2^ 



(n + 2)(n + /3 £ + 3)(n + 2/^ + 3) 1 



(n + A + 2) 



2^ 



(7l-l)(7l + ^-l)(7l + 2^ ) 

(n + pi) 



where 



4fiD r 



h 2 (n + /3 e + l) 



(54) 
(55) 



Table 1: Reduced masses and spectroscopic properties of the CO and NO diatomic 
molecules in the ground electronic state. The data listed in this table are taken from 

[23 ES| . 



Parameters 


CO 


iVO 


A> (in eV) 
r G (in ,4°) 
//(in amu) 


10.84514471 

1.1282 
6.860586000 


8.043782568 

1.1508 
7.468441000 
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Table 2: Comparison of the energy eigenvalues (in eV), corresponding to the Kratzer 
potential for various n and I quantum numbers for CO and NO diatomic molcules, 
where he = 1973.29elM°. 







CO [FM] 


CO [EQR] 


CO [AIM] 


NO[FM] 


JVO[EQR] 


NO [AIM] 


o 


o 


-10 


79431534387622 


-10 


794315323 


-10 


79431532 


-8 002658755212952 


-8 


002659419493 


-8 


002659417 


1 





-10 


69383913769446 


-10 


693839925 


-10 


69383992 


-7.921456003136883 


-7 


921456840689 


-7 


921456839 


1 


1 


-10 


69337109882925 


-10 


693371229 


-10 


69337123 


-7.921042972272428 


-7 


921043829925 


-7 


921043834 


2 





-10 


59476059512928 


-10 


594760890 


-10 


59476089 


-7.841483226715529 


-7 


841483958093 


-7 


841483956 


2 


1 


-10 


59429734443464 


-10 


594298692 


-10 


59429869 


-7.841075958492119 


-7 


841077185904 


-7 


841077188 


2 


2 


-10 


59337288634942 


-10 


593374417 


-10 


59337441 


-7.840262393682968 


-7 


840263768523 


-7 


840263771 


3 





-10 


49705105930509 


-10 


497052462 


-10 


49705246 


-7.76271486692960 


-7 


762716067159 


-7 


762716066 


3 


1 


-10 


49659537021400 


-10 


496596643 


-10 


49659664 


-7.762314093207729 


-7 


762315408528 


-7 


762315413 


3 


2 


-10 


49568316746262 


-10 


495685124 


-10 


49568512 


-7.761512650377457 


-7 


761514215884 


-7 


761514218 


3 


3 


-10 


49431563207918 


-10 


494318144 


-10 


49431814 


-7.760311623981464 


-7 


760312738370 


-7 


760312744 


4 





-10 


40068823210378 


-10 


400689478 


-10 


40068947 


-7.685127711083608 


-7 


685129080626 


-7 


685129079 


4 


1 


-10 


40023809235447 


-10 


400239921 


-10 


40023992 


-7.684732464238659 


-7 


684734413653 


-7 


684734417 


4 


2 


-10 


39933885575065 


-10 


399340924 


-10 


39934092 


-7.683942903888297 


-7 


683945202003 


-7 


683945203 


4 


3 


-10 


39798982948594 


-10 


397992722 


-10 


39799272 


-7.682760084230797 


-7 


682761690175 


-7 


682761696 


4 


4 


-10 


39619321456787 


-10 


396195666 


-10 


39619567 


-7.681181936089156 


-7 


681184244677 


-7 


681184246 


5 





-10 


30564563709163 


-10 


305647347 


-10 


30564735 


-7.608697429711712 


-7 


608699510108 


-7 


608699509 


5 


1 


-10 


30520186777971 


-10 


305203938 


-10 


30520394 


-7.608308414224594 


-7 


608310715917 


-7 


608310719 


5 


2 


-10 


30431444380035 


-10 


304317236 


-10 


30431723 


-7.607531297904489 


-7 


607533247563 


-7 


607533248 


5 


3 


-10 


30298450330134 


-10 


302987469 


-10 


30298747 


-7.606365523492213 


-7 


606367345012 


-7 


606367349 


5 


4 


-10 


30121238898166 


-10 


301214985 


-10 


30121499 


-7.604810653591787 


-7 


604813367976 


-7 


604813368 


5 




-10 


29899764904374 


-10 


299000242 


-10 


29900024 


-7.602869548627632 




602871795644 




602871795 
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Table 3: Comparison of the energy eigenvalues (in eV), corresponding to the modi- 
fied Kratzer potential for various n and I quantum numbers for CO and NO diatomic 
molcules, where he = 1973.29eVM°. 



n 




U(J [r 1VLJ 


L (_/ [rji^rCJ 


O Kj [IN U J 


i V (J [r 1V1J 


iv L7[rjV^KJ 


AT/~I r\TT Tl 
iV C[1N UJ 








0.050828927974365U0 


0.050829386733 


0.050823 


0.04112347897894253 


0.041123148507 


0.041118 


1 





0.1513051341561269 


0.151304784801 


0.151287 


0.1223262310550117 


0.122325727312 


0.122311 


1 


1 


0.1517731730213381 


0.151773481462 


0.151755 


0.1227392619194667 


0.122738738076 


0.122724 


2 





0.2503836767213095 


0.250383819984 


0.250354 


0.2022990074763653 


0.202298609907 


0.202274 


2 


1 


0.2508469274159442 


0.250846018075 


0.250816 


0.2027062756997760 


0.202705382097 


0.202681 


2 


2 


0.2517713855011650 


0.251770292931 


0.251744 


0.2035198405089265 


0.203518799478 


0.203494 


3 





0.3480932125454999 


0.348092247640 


0.348051 


0.2810673574989346 


0.281066500842 


0.281033 


3 


1 


0.3485489016365868 


0.348548066746 


0.348507 


0.2814681409841651 


0.281467159473 


0.281434 


3 


2 


0.3494611043879683 


0.349459585720 


0.349418 


0.2822695838144371 


0.282268352117 


0.282235 


3 


3 


0.3508286397714091 


0.350826566166 


0.350785 


0.2834706102104301 


0.283469829631 


0.283436 


4 





0.4444560397468020 


0.444455232045 


0.444403 


0.3586545231082869 


0.358653487375 


0.358611 


4 


1 


0.4449061794961153 


0.444904789014 


0.444852 


0.3590497699532351 


0.359048154348 


0.359006 


4 


2 


0.4458054160999403 


0.445803785755 


0.445751 


0.3598393303035978 


0.359837365998 


0.359795 


4 


3 


0.4471544423646421 


0.447151987956 


0.447099 


0.3610221499610970 


0.361020877826 


0.360978 


4 


4 


0.4489510572827147 


0.448949044337 


0.448895 


0.3626002981027385 


0.362598323324 


0.362555 


5 





0.5394986347589601 


0.539497362596 


0.539434 


0.4350848044801827 


0.435083057893 


0.435032 


5 


1 


0.5399424040708798 


0.539940771611 


0.539877 


0.4354738199673003 


0.435471852084 


0.435421 


5 


2 


0.5408298280502333 


0.540827474443 


0.540764 


0.4362509362874052 


0.436249320438 


0.436198 


5 


3 


0.5421597685492472 


0.542157240772 


0.542093 


0.4374167106996811 


0.437415222989 


0.437364 


5 


4 


0.5439318828689306 


0.543929725307 


0.543865 


0.4389715806001080 


0.438969200025 


0.438917 


5 


5 


0.5461466228068463 


0.546144468004 


0.546082 


0.4409126855642622 


0.440910772357 


0.440858 
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With equations ( |53l) and ( )54l) . we can deduce the following relations: 



Tx = i{R n/ (r)\r\ R n/ (r)) + (r) | r— | (r)) 



(n + ft ) - 



2\ 



(n + 2)(n + ft + 3)(n + 2ft + 3) _ 3 
(n + ft + 2) 2\ 



(n-l)(n + ft-l)(n + 2ft ) 



(n + ft ) 



(56) 



and 



T 2 = t{R n ,t (r)\r\ R n/ (r)> - (i^ (r) | r— | (r)) 



(n + ft +1) 



2\ 



(n + 2)(n + ft + 3)(re + 2ft + 3) _ 3 
(n + ft + 2) 2\ 



(w-l)(7i + ft-l)(ra + 2ft ) 
(n + ft ) 



(57) 

In this case, the matrix elements for the two potentials (Kratzer and modified Kratzer 
potentials) give the same results, since these matrix elements depend on b = DqTq and 
a = —2D Q r only, and not on c = D Q . Hence, when the Kratzer potential is shifted by 
the amount D , (that is, the modified Krazter potential), it has no effect on the matrix 
elements. For the numerical results, see Tables 4 and 5. 
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Table 4: The values of the (R n> e \ r \ R Hj e) [Equation (53)], (R n> £ \ r-^ \ R Hj 
[Equation (54)], Ti [Equation (56)] and T 2 [Equation (57)], corresponding to the Kratzer 
and modified Kratzer potentials for various n and £ quantum numbers for CO diatomic 
molcule, where he = 1973.29elM°. 



n 




(-Rn, f I r | ii n , ^> 


(R n , e\r£\ Rn, l) 


Ti 


r 2 


1 





-0 


4761490924054464 


16.97323506451677 


195 


3895273866551 


161 


4430572576216 


1 


1 


-0 


4761710273791221 


16.97343118040934 


195 


3940309777938 


161 


4471686169751 


2 





-0 


4107558210705384 


9.473123928272610 


162 


6714457007481 


143 


7251978442029 


2 


1 


-0 


4107755614421622 


9.473236107448747 


162 


6755771602085 


143 


7291049453111 


2 


2 


-0 


4108150436563324 


9.473460462197114 


162 


6838400978704 


143 


7369191734761 


3 





-0 


3789129917968390 


7.666855095230551 


148 


3356842660929 


133 


0019740756318 


3 


1 


-0 


3789316534237137 


7.666946649417545 


148 


3396482090482 


133 


0057549102131 


3 


2 


-0 


3789689781099853 


7.667129754891770 


148 


3475761191884 


133 


0133166094049 


3 


3 


n 
-u 




7.667404405854176 


148 




-LOO 




4 





-0 


3533372017343955 


6.615793772861618 


137 


1863080576614 


123 


9547205119381 


4 


1 


-0 


3533549972979473 


6.615873160774083 


137 


1901403213880 


123 


9583939998398 


4 


2 


-0 


3533905898264022 


6.616031934119029 


137 


1978048773512 


123 


9657410091131 


4 


3 


-0 


3534439821224761 


6.616270087936851 


137 


2093017825661 


123 


9767616066924 


4 


4 


-0 


3535151205980674 


6.616587357061658 


137 


2246186798360 


123 


9914439657127 


5 





-0 


3312181486228352 


5.899103372782779 


127 


7353299114129 


115 


9371231658474 


5 


1 


-0 


3312351970446675 


5.899174359294310 


127 


7390497423004 


115 


9407010237118 


5 


2 


-0 


3312692952626132 


5.899316330129842 


127 


7464894362121 


115 


9478567759524 


5 


3 


-0 


3313204460252490 


5.899529280914656 


127 


7576490574162 


115 


9585904955868 


5 


4 


-0 


3313885980893305 


5.899812974634305 


127 


7725166245981 


115 


9728906753295 


5 


5 


-0 


3314738123034079 


5.900167633018892 


127 


7911043436078 


115 


9907690775700 
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Table 5: The values of the (R n> e \ r \ R Hj e) [Equation (53)], (R n> £ \ r-^ \ R Hj 
[Equation (54)], Ti [Equation (56)] and T 2 [Equation (57)], corresponding to the Kratzer 
and modified Kratzer potentials for various n and £ quantum numbers for NO diatomic 
molcule, where he = 1973.29elM°. 



n 


I 


(_R nj ^ | r | ii nj ^) 


(R n , t\r£\ Rn, l) 


Ti 


r 2 


1 





-0.4819553974463722 


16.21611432627930 


178.3853413621973 


145.9531127096387 


1 


1 


-0.4819818755489971 


16.21633768030118 


178.3902449613004 


145.9575696006980 


2 





-0.4124035377598499 


9.040224937878095 


147.1042625831576 


129.0238127074014 


2 


1 


-0.4124272376527343 


9.040352590403065 


147.1087421492765 


129.0280369684704 


2 


2 


-0.4124746395526131 


9.040607890585319 


147.1177013069172 


129.0364855257466 


3 





-0.3785773170199991 


7.313655510582683 


133.4152241892145 


118.7879131680491 


3 


1 


-0.3785996534470747 


7.313759628367830 


133.4195129653269 


118.7919937085912 


3 


2 


-0.3786443283574571 


7.313967860029029 


133.4280905503876 


118.8001548303296 


3 


3 


-0.3787112793749515 


7.314279887895779 


133.4409442456685 


118.8123844698770 


4 





-0.3514058390444385 


6.309618628486895 


122.7741267465619 


110.1548894895881 


4 


1 


-0.3514270813879861 


6.309708856787779 


122.7782656025786 


110.1588478890030 


4 


2 


-0.3514695680843574 


6.309889310052519 


122.7865433532380 


110.1667647331330 


4 


3 


-0.3515332399199089 


6.310159713089281 


122.7989477575504 


110.1786283313719 


4 


4 


-0.3516182293790348 


6.310520592939504 


122.8155035676796 


110.1944623818006 


5 





-0.3278978812775772 


5.625407532078025 


113.7570644421939 


102.5062493780379 


5 


1 


-0.3279181800627655 


5.625488165029850 


113.7610753342060 


102.5100990041463 


5 


2 


-0.3279587796016631 


5.625649427995388 


113.7690971617573 


102.5177983057665 


5 


3 


-0.3280196234075503 


5.625891075135918 


113.7811180744773 


102.5293359242055 


5 


4 


-0.3281008382258643 


5.626213577576859 


113.7971620775729 


102.5447349224192 


5 


5 


-0.3282023110716663 


5.626616443662535 


113.8172054699764 


102.5639725826513 



5 Conclusions 

In this paper, I have studied the eigenvalues, the eigenfunctions and the matrix elements 
of the Kratzer- type molecular potentials (generalized Kratzer molecular potentials). The 
ladder (creation and annihilation) operators for the radial wavefunctions are established. 
Also, the Hermitian operators of these ladder operators are obtained. These operators 
satisfy the commutation relations of an S77(l,l) dynamical group, the action of £± 
on the wavefunctions that reveals the fact about the raising and lowering effect are 
established. The matrix elements of the different functions r and r4- are also obtained 

dr 

from the ladder operators in a closed form. 

The solutions of the generalized Kratzer molecular potentials are obtained via an 
SU(1, 1) algebraic approach and the results can be generalized to other form of potentials 
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in equation (jSJ). The generalized Kratzer potential model proposed in this work [T8| IT9]. 
allows one to obtain the eigenvalues, the matrix elements and the radial eigenfunctions 
for the two important molecular potentials simultaneously (that is, the Kratzer and 
modified Kratzer potentials). 

This generalization covers the descriptions about the two potentials, the shapes of 
the two potentials for CO and NO diatomic molecules are also shown in figure 1. More 
importantly, the generalized Kratzer potential model used in this work generates results 
for the two molecular potentials. In addition, the numerical results are obtained for the 
eigenvalues for CO and NO diatomic molecules for the two molecular potentials. The 
results obtained are compared with other existing results (EQR, AIM and NU), and the 
results agree favourably with other results, see Tables 2 and 3. 

Tables 4 and 5 show the matrix elements of r, r4-, their sum (1^) and their difference 
(T2) obtained for CO and NO diatomic molecules, the results obtained are the same 
results for the two molecular potentials. This is because, the matrix elements depend 
on b = DqTq and a = — 2Do r o only, and not on c = Dq. Hence, the shifted amount Dq, 
(in the case of the modified Krazter potential) has no effect on the matrix elements. 

The advantage of the present approach is that it enables one to find the energy eigen- 
values, eigenfunctions and matrix elements in a simple way. The approach presented in 
this study is efficient and is a very useful link for finding the matrix elements from lad- 
der operators. This approach can be used to find the energy eigenvalues, eigenfunctions 
and radial matrix elements of the Schrodinger equation with a given exactly solvable 
molecular potentials for various diatomic molecules for any values of n and i quantum 
numbers. 

The ladder operators constructed in this study, are very useful tools in quantum- 
mechanical calculations of the various matrix elements based on the Kratzer molecular 
basis function. Furthermore, these operators can be used in constructing coherent states. 

[Note to the desk Editor: Please, put the figure in the section 2, where it is indicated 
in the text] 
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